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Abstract 

We present an extensive but concise review of our present understanding, largely 
based on theory and simulation work from our group, on the equilibrium behavior 
of solid surfaces and nanosystems close to the bulk melting point. In the first part 
we define phenomena, in particular surface melting and nonmelting, and review 
some related theoretical approaches, from heuristic theories to computer simula- 
tion. In the second part we describe the surface melting/nonmelting behavior of 
several different classes of solids, ranging from van der Waals crystals, to valence 
semiconductors, to ionic crystals and metals. In the third part, we address special 
cases such as strained solids, the defreezing of glass surfaces, and rotational sur- 
face melting. Next, we digress briefly to surface layering of a liquid metal, possibly 
leading to solid-like or hexatic two dimensional phases floating on the liquid. In the 
final part, the relationship of surface melting to the premelting of nanoclusters and 
nanowires is reviewed. 
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1 Introduction: surface melting/nonmelting 

Melting, certainly one of the longest known phase transitions, always appeared 
a very special one. Although undoubtedly first order, it curiously displays 
only half the hysteresis cycle one would expect. In a regular first order phase 
transition, the free energies of two phases cross, while both remain locally 
stable beyond the crossing. As a result, it is generally possible to undercool the 
high temperature phase, and to overheat the low temperature phase. Periodic 
heating and cooling through a first order transition will thus generate a full 
hysteresis cycle. 

In melting, this is only half true. It is generally possible and easy to supercool a 
liquid: a glass of water may not freeze at all during an icy night, only to do that 
suddenly when disturbed the next morning. On the contrary, it is surprisingly 
difficult and often impossible to overheat a free solid. The reason for that is 
that a free solid has surfaces; and when the melting temperature is approached, 
melting of the bulk crystal is ready to begin from there. Supercooling of the 
liquid is allowed by the absence of a solid germ; overheating of the solid is 
prevented by that omnipresent germ, a wet surface. This observation provides 
perhaps the first macroscopic hint that solid surfaces might be wet already 
somewhat below the melting point. The starting motivation for a microscopic 
study of surface melting begins right there. 

Consider a semi-infinite homogeneous solid, e.g., a crystal with a free surface, 
in thermodynamic equilibrium at temperature T and pressure P with its own 
vapor (Fig. 1). As T and P are raised to approach melting at the bulk triple 
point (T m , P m ), where solid, gas and liquid coexist, the solid- vapor interface 
will generally wet itself with an atomically thin liquid film, whose thickness 
£(T) — > oo when T -> T m . This is called surface melting (SM). A rather 
famous microscopic characterization of SM is that of Pb(llO), first reported 
by Frenken and van der Veen [1] . Medium energy ion scattering data reveal 
(Fig. 2) the presence at the solid- vapor interface of a disordered (liquid) film 
whose thickness £(T) grows without limit as T — > T m . A number of case 
studies are now known in the literature, and SM turns out to be the common 
behavior for many surfaces. It is however not general, and many surfaces are 
known that do not melt in the same sense, but remain solid and crystalline 
all the way up to T m . Not surprisingly, this behavior has been called surface 
nonmelting (NM) [2]. Slightly more surprisingly, SM and NM can both occur 
for the same substance, of course on two different crystallographic faces. For 
example, NM was demonstrated for Pb(lll) [3], while by contrast SM prevails 
on Pb(llO) [1]. 

In the language of wetting, SM is nothing else than triple point wetting[4], 
namely complete wetting of the solid-vapor interface by the liquid, in other 
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words wetting of the solid by its own melt. By contrast, NM corresponds to 
partial wetting, schematized in Fig. 3. In this case, there is a clear connec- 
tion between the surface thermodynamic parameters and the angles of partial 
wetting [5] as will be clarified further below. Thermo dynamically, NM will 
take place whenever the solid surface (solid- vapor interface) is unable to lower 
its free energy by turning into a sequence of two separate solid-liquid plus 
liquid-vapor interfaces, namely when 



7sv<7sl + 7lv, (1) 



where the 7's denote the free energies of the three solid- vapor (SV), solid-liquid 
(SL), and liquid- vapor (LV) interfaces [7]. Conversely, when the inequality (1) 
does not hold, SM will ensue instead of NM. As Fig. 2 shows, in SM the liquid 
film, although technically of divergent thickness at T m , remains atomically thin 
until fractions of a degree below T m . Interactions propagate readily across such 
a thin liquid film. The solid-liquid and the liquid-vapor interfaces "feel" their 
mutual presence, so much that they actually merge when £{T) is small enough. 
For a generic temperature T < T m the free energy change caused by melting 
a solid film into a liquid one of thickness £ can be written 



AG I A = [( 7S L + 7lv) - 7sv] + £L Pl (T m -T)+ V{£) (2) 



where A is the surface area, L is the latent heat of melting per unit mass, p\ is 
the density of the liquid, and V(£) represents a phenomenological "interaction" 
free energy term between the liquid-vapor and the solid-liquid interfaces. As 
defined, the interaction V(£) vanishes for £ — > 00, and tends to V(0) = 7sv — 
(7sl + 7lv) when £ — > 0, when the dry solid surface free energy is recovered. 
Eq. (1) amounts to say that in SM the interface interaction V(£) is globally 
repulsive, whereas in NM it is attractive. A useful mnemonic is that SL-LV 
repulsion causes the liquid film to expand (SM), while SL-LV attraction causes 
the liquid film to collapse and disappear (NM). 

The microscopic mechanism behind interface attraction or repulsion varies, 
and may be eventually related to such different physical causes as layering 
order [6] , van der Waals forces [8] , or molecular order [9] . We shall provide 
specific cases further below. 
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2 Theory; brief review 

2.1 Heuristics: solid surface instability 

Early crude but appealing reasoning leading to SM has been based on heuris- 
tics. Bulk melting usually correlates well with the so called Lindemann melt- 
ing criterion: when the r.m.s. thermal vibration amplitude of the solid reaches 
some 14% of the nearest-neighbor distance [10], solids generally melt. Atoms 
at surfaces are less coordinated and generally shakier than in bulk, and their 
vibration amplitudes are correspondingly larger. Surfaces will thus reach the 
Lindemann instability vibration amplitude at a lower temperature than the 
bulk. A model embodying the bulk Lindemann criterion is the mechanical 
thermal instability model [11]. An infinite solid will, if prevented from melting, 
become mechanically unstable at a sufficiently high temperature. Although the 
ideal mechanical instability temperature of a solid is different, and of course 
somewhat higher than the true melting temperature T m (where the free energy 
crossing of solid and liquid phases takes place), it can nonetheless heuristically 
be taken as a qualitative indicator of the tendency of the solid to melt. Adopt- 
ing that line, one can ask whether the stability of the surface will not cease 
before that of the bulk. In a model semi-infinite solid it is found indeed that the 
mechanical instability of the first surface layer occurs at a temperature which is 
only 73% that of the bulk [12]. This suggests - correctly - that surfaces should 
begin to soften and melt at about 3/4 the bulk melting temperature. That 
had been long known and was noted by Tammann and Stranski [13,14,15]. 
Self-consistent surface phonon calculations [16] and also experiments [17] later 
revealed a very pronounced anharmonic outward expansion of the first surface 
layer relative to that of the bulk. Although it is not really possible to gener- 
alize, this is undoubtedly one of the qualitative elements heralding a stronger 
tendency of a surface to become unstable before the bulk. 



2.2 Density functional theory of surface melting 

The simplest proper theory of SM is based on optimizing the order parameter 
profile, in particular the density profile p(x) of a solid-vapor interface in a sys- 
tem with short range interactions. Assuming a bulk grand potential density 
u = u(p), the solid vapor coexistence near the triple point as in Fig. 1 implies 
two equivalent solid and vapor minima of the grand potential per unit vol- 
ume uj{p) plus one additional secondary liquid minimum with a higher grand 
potential density. The total grand potential of the inhomogeneous system can 
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be written 

n = JJJ d3x 



c( P (x)) + J - \v P y 



(3) 



In this Ginzburg-Landau (or Cahn-Hilliard) type phenomenological free en- 
ergy form (acceptable for a system with short range forces) the first term 
represents the thermodynamic potential density for a uniform system, while 
the gradient term signifies the extra cost caused by any spatial change in the 
order parameter. For a solid-vapor interface [18] we may assume p = p{z) and 
seek to minimize Q (per unit area): 



+00 

n/A = J c 



dz 



u(p(z)) + 




mm 



(4) 



with the constraints p(— 00) = p s and p(+oo) = p v , where p s and p v are 
solid and vapor densities respectively. Amusingly, this minimization problem 
is formally identical to the minimum action problem for the one dimensional 
lagrangian motion of a classical point object of coordinate p, mass J, as a 
function of time z, in a potential — uj{p). Starting at time z = —00 with zero 
kinetic energy and with a small p = p v from a first hilltop of height — cu(p v ) 
(the vapor phase), the point moves "downhill" eventually reaching at time 
z = +00 the last, and exactly equivalent, hilltop — oo(p s ) = —u(p v ) (the solid 
phase) (Fig. 4). 

At temperatures well below the melting point, there is no liquid phase (not 
even metastable), meaning that there exists no liquid-like local minimum in 
ui{p). Thus — u(p) has no other maxima than the vapor and solid hilltops, 
and the resulting solid-vapor interface is unsplit and featureless. Close to T m 
however — u(p v ) develops a secondary maximum — a lower hilltop, as it were— 
at an intermediate p = p\. En route from the vapor to the solid hilltops, 
the point particle must negotiate this intermediate "mountain pass" , where it 
will not stop, but will still slow down considerably. It is easy to check that if 
A = — cu(ps) — (-ui(pi)) is the height difference between main and secondary 
hilltops (the free energy difference between solid and liquid, that goes to zero 
only at the melting point), then the time 5z the point will spend near the liquid 
hill, and therefore with p w pi is proportional to log A. Thus the interface 
density profile is now split into two interfaces: first a vapor-liquid one, then 
a liquid-solid one. They are separated by a liquid film of density close to pi, 
whose thickness £ = 5z is logarithmically increasing as T — > T m . And this is 
of course just surface melting. 

The above also suggests that one might identify some temperature T w < T m 
where the bulk free energy first develops the local minimum corresponding 
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to the liquid phase, with the temperature T w where a solid surface is likely 
to begin wetting itself with the thinnest liquid film. It must be underscored 
that this simple theory only describes SM and cannot account for NM or for 
many other complications, upon which we shall return below. The presence of 
a "surface term" in Eq.(4) will additionally modify this scenario [18]. 

2.3 Lattice mean field theory of surface melting (Trayanov-Tosatti), and be- 
yond 

The discussion above is centered on density as the sole order parameter. A 
solid actually differs from the liquid not just by the average atomic density 
p , but also by the infinite set of crystalline Fourier components of the atomic 
density pc (G denoting reciprocal lattice vectors) that are nonzero only in the 
solid. A better order parameter theory of SM than that outlined above should 
be based on a free energy expression that correctly includes these crystalline 
order parameters. Trayanov and Tosatti [19] built a simplified but microscopic 
two-order parameter lattice theory based on the average density p P m s a sec- 
ond "crystallinity" order parameter c, the latter effectively replacing the less 
manageable infinite set of pc- Through a further mean field approximation 
where lattice layers are assumed to behave uniformly (thus neglecting rough- 
ening fluctuations), realistic systems such as Lennard- Jones (LJ) solid surfaces 
can be described in this manner. Fig. 5 illustrates the surface melting of fee 
LJ(llO) described by this lattice theory. The approach also demonstrated how 
in the LJ systems the logarithmic liquid film thickness divergence close to T m 
is asymptotically replaced by a more realistic (T m — T) -1 / 3 divergence, due to 
the (omnipresent) van der Waals long range interatomic potential tails. Let 
us expand briefly on this point. Given three media S, L and V, with a finite 
thickness £ of L sandwiched between semi-infinite S and V, there will generally 
arise a long range tail to the interaction free energy between the two interfaces 

V(l) ~ H/e (5) 

due to long range dispersion forces. The so-called Hamaker constant H is 
generally positive when the intermediate L phase is dielectrically less dense 
that S, and more dense than V[20]. Positive H means long range repulsion 
between the interfaces, which favors SM. That is precisely the case for the 
LJ-like systems, where the liquid is less dense than the solid. The same will 
generally also apply to many metals and other materials. In different cases, 
such as valence semiconductors and semimetals, melting leads to a denser 
and generally metallic liquid state, and there H < 0. Negative expansion at 
melting and a negative Hamaker constant occur also in other notable cases, 
water among them [21,22]. Here H < means a long range attraction between 
interfaces, and that implies NM. As we shall see below, the surfaces of a 
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semiconductor actually appear to wet close to T m . However the wetting is 
incomplete, and the melted film thickness £ remains microscopically thin, in 
agreement with their negative Hamaker constant. 

Of course mean field theories, such as Trayanov-Tosatti above, are only an 
approximation to a real critical phenomenon. Mean field is usually a bad 
approximation sufficiently close to a critical point (in SM the critical point is 
T = T m ), where the correlation length (here signified by the liquid thickness 
t) diverges. In systems with short range forces critical fluctuations may lead 
to nonclassical exponents. For SM, this aspect was discussed theoretically 
by Lipowsky and collaborators [18], and later by Chernov and Mikheev[23]. 
It was also addressed through atomistic simulations by Chen et al [24] . The 
conclusion is that although non-mean field fluctuations (mostly roughening 
fluctuations) actually have a large quantitative effect [24], they will not change 
the mean field exponent —1/3 which originates from long range forces [8]. 

2.4 Simulations: Molecular Dynamics 

Atomistic simulations constitute a very important tool in the field of melting 
and surface melting. They can be either of Monte Carlo (MC) type, or of 
Molecular Dynamics (MD) type. We shall concentrate on the latter, using 
work done in our group for specific examples. SM simulations are generally 
conducted on crystal slabs, made up of a sufficiently large number of atomic 
solid layers, with two free surfaces, or with one free and one "frozen" (fixed 
atoms) surface, and periodic boundary conditions along the two directions 
parallel to the surfaces. The in-plane simulation cell size must be allowed or 
otherwise made to expand gradually as temperature increases, to guarantee 
that the planar stress in the solid on account of thermal expansion remains as 
close as possible to zero. The simulation consists of solving Newton's equations 
of motion 

= -Vi-E P ot({R}) (6) 

for all particles of mass Mj, and coordinates Rj. The basic input needed is 
the potential energy E pot ({H}) as a function of all coordinates. Depending 
on systems, and on the accuracy needed, the choice is between empirical in- 
teratomic potentials on one side, and first principles total energy calculations 
£ , pot ({R}) on the other side. Empirical two-body potentials, such as Lennard- 
Jones (LJ), are quite reasonable for rare gas solids or van der Waals molecular 
crystals, but not for most other solids. For metals, empirical many body po- 
tentials such as the Embedded Atom Model [25], the similar glue model [26], 
the Finnis-Sinclair potential [27] etc., are much more suitable than two-body 
forces. Many-body empirical potentials were developed also for such systems 
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as valence semiconductors [28,29] but here the first principles simulation ap- 
proach is generally much more appropriate [30] . Typical SM simulation density 
profile outputs are shown in Fig. 6 for a modified LJ system and for a Au metal 
surface respectively. Technical details and some further results for simulations 
of SM can be found in the review article by Di Tolla et al [31]. Simulations 
may in many ways replace experiment, and must be similarly regarded and 
understood by means of theory — not discounted as self-explanatory as it is 
sometimes tempting to do. Similar to experiment, they do uncover (occasion- 
ally by serendipity), novel or unexpected behavior. Such has been the case for 
surface nonmelting, first discovered in simulation [2] and next found experi- 
mentally [3] , and several other phenomena, to be discussed further below. 



3 Surface melting and related phenomena in real solids 

3. 1 van der Waals solids: interplay of roughening, preroughening , and melt- 
ing 

Rare gas solids are reasonably described through a Lennard- Jones potential. 
As first shown by extensive MD studies by Broughton and collaborators [32,33] 
they are expected to display surface melting on all crystalline faces. Here truly 
microscopic experiments are not abundant, because of technical difficulties. 
There is however a large body of work by the surface adsorption community, 
interested in understanding adsorption isotherms as a function of tempera- 
ture [34]. 

Adsorbtion of atoms on a substrate generally takes place in a layer-by-layer 
mode, up to the so-called surface roughening temperature Tr, where the inter- 
face width diverges and layers cease to exist. Surface roughening-a Kosterlitz- 
Thouless phase transition taking place (sometimes but not always) at a tem- 
perature Tr < T m — is a totally distinguished phenomenon from surface melt- 
ing [35], but it does have some relations with it. It can be argued for example 
that surfaces that undergo SM at T = T m will necessarily undergo a roughen- 
ing transition at T R < T m [36]. That is because a liquid surface is technically 
rough, so that capillary fluctuations cause the liquid-vapor interface width to 
diverge. Conversely, although there is apparently no theorem that says so, a 
NM surface will usually not undergo roughening below T m . Disregarding con- 
nections with surface melting — even though they are sometimes important, as 
discussed below — roughening is generally discussed in purely lattice models. 

There is yet another surface phase transition that may occur below melting, in 
fact below roughening, therefore called preroughening (PR) [37,38,39,40]. PR 
is associated with the bulk lattice planes having some stacking such as abab . . . 
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or abcabc . . . where planes a and b can only be made to coincide through a frac- 
tional translation. At low temperatures a perfect surface, say a (100) face, will 
with identical probability be either a or b terminated. As in an Ising model, 
either choice will "spontaneously break" the symmetry between a and b . At 
the PR temperature Tpp, a surface restores precisely the a-b symmetry in an 
Ising-like fashion. To achieve that, the topmost atomic layer undergoes a most 
surprising spontaneous rearrangement, from full coverage below Tpp, (full a 
or full b) to | coverage above T PR . The excess atoms simply migrate away, 
turning the top half layer into a sort of soft a-b "checkerboard". Strictly at 
T = T PR the surface becomes technically rough, only to become flat again be- 
tween Tpp and Tp. Preroughening was first discovered theoretically in lattice 
models [37], and received little attention by experimentalists until the equally 
surprising phenomenon of re-entrant layering was discovered [41,42]. Adsorp- 
tion of Ar, Xe, Kr on a graphite substrate (Fig. 7) proceeded by full layers 
below 0.82 T m . This appeared to be roughening: but layering surprisingly re- 
entered at higher temperature, this time by half layers. While that clearly 
called for an explanation in terms of preroughening [43] a considerable con- 
troversy arose, as independent evidence suggested that the first surface layer 
was in reality no longer solid, with considerable liquid-like disorder and first 
layer mobility in that temperature regime [44] . Moreover, reentrant layering 
signaled a clear first order transition, whereas PR should be second order. 
The resolution of this puzzle, reached through theory [45], grand canonical 
MC [46,47] and MD [48] simulations, proved to be that PR is linked to first 
(half-) layer surface melting: in essence both phenomena take place simulta- 
neously. PR in a rigid lattice surface model would indeed only take place at 
higher temperatures, and the half monolayer would in that case be stabilized 
through a second order transition. In real off-lattice Argon however the half 
monolayer soft checkerboard immediately melts, with a first order transition 
and a large positional entropy gain over the corresponding hypothetical solid 
half monolayer. As a result, preroughening and first layer melting actually 
take place simultaneously, at a much lower temperature than either of them 
separately, and of course through a first order transition. Fig. 8 shows a realis- 
tic snapshot picture, obtained by MC, of the Ar(lll) surface simulated in the 
reentrant layering temperature regime, showing the spontaneous formation of 
a roughly half-coverage, liquid like top monolayer. 



3.2 Ionic insulators 



Surfaces of ionic crystals such as NaCl, other alkali halides, MgO, etc, are 
often used as substrates for growing other materials, and their properties are 
therefore better known far below the melting point. Yet, some data are also 
available about their behavior at T m . 
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In this respect, NaCl appears to be one of the best studied, and can therefore 
be chosen as a case study. Argon bubble studies of liquid NaCl in contact with 
the solid revealed a surprising lack of complete wetting, with a large partial 
wetting angle of about 48 degrees (Fig. 9) [49] . Very recent MD simulations of 
solid NaCl(lOO) clearly demonstrates NM of this solid surface, also predicted 
to survive in a metastable state well above T m [9]. Simulations of a droplet of 
melted NaCl brought into contact with solid NaCl(lOO) at T = T m (Fig. 10) 
clearly shows incomplete wetting, with an external partial wetting angle of 
(50 ± 5)°, in good agreement with the bubble experiments [50]. 

The microscopic reasons for the exceptional stability of the solid NaCl(lOO) 
surface and its unreadiness to wet itself are presently under scrutiny. They are 
intriguing in view of a large positive Hamaker constant - NaCl undergoing a 
26% volume expansion at melting - which would suggest the opposite behav- 
ior, namely SM. A second element wrongly suggesting SM against NM is the 
structure of the liquid surface, very strongly fluctuating and totally devoid of 
layering (Fig. 11) unlike e.g., nonmelting metal surfaces (see below). So why 
does strong NM arise in alkali halides? 

Preliminary useful clues come from surface thermodynamics. MD simulations 
permit the explicit calculation of ^lv-iIsl and '-fsv as a function of T[9]. It 
appears that while '-fsv drops very strongly in virtue of large surface relax- 
ations and the associated growing vibrational entropy at T ~ T m , the liquid 
surface is incapable of developing a competing amount of entropy, presumably 
because of the all-important local charge neutrality constraint, and of the re- 
lated appearance of molecular order in the fluctuating liquid NaCl surface. 
This physical result is the subject of current work [96]. 



3.3 Valence semiconductors and semimetals 

Valence semiconductors such as Si, Ge, GaAs and semimetals like Ga are 
known to turn fully metallic and denser either at high pressure or when they 
melt at high temperature. This is illustrated by the schematic Ga phase dia- 
gram of Fig. 12. The higher liquid density here causes the Hamaker constant 
to be negative [24] , and this necessarily hinders complete surface melting. Re- 
ports of regular SM with unlimited growth of the liquid film at T m that have 
appeared for Ga [51] are incompatible with the negative Hamaker constant, 
and might in our view be an artifact, possibly due to strains (see below). Due 
to the long range attraction, regular and complete SM is excluded for all faces 
of these H < materials. 

The clean surfaces of semiconductors, Ge(lll) and Si(lll), were thoroughly 
investigated for surface melting. While they were found to disorder [52] and 
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simultaneously to become metallic [53] on the topmost layer at some T < T m , 
there is clear evidence that the disordered film does not grow as T m is ap- 
proached. In order to understand theoretically the top layer disordering and 
metallization, MD simulations are perhaps the best tool. Since electronic ef- 
fects and the semiconductor-metal transition play a very important role in 
the melting of these systems, it seems mandatory to take recourse in this case 
to first-principles simulations, where forces acting on the ions are obtained 
by solving the full electronic structure problem. Simulations of this type, pi- 
oneered by Car and Parrinello in the mid 80s [54] , are by their nature much 
more expensive and demanding than the empirical ones. Results for Ge(lll) 
by Takeuchi et al [30] are shown in Fig. 13 at T ~ T m . As signaled by the 
broader density profile, the topmost bilayer is indeed strongly disordered, but 
disorder does not propagate to the second bilayer, indicating incomplete melt- 
ing. Electronically, the surface disordered film is metallic [53] , as opposed to 
the underlying bulk, which remains semiconducting. 



3.4 Metals 



The surface melting behavior of metallic surfaces has received by comparison 
much more experimental attention. Metal surfaces usually possess generally 
positive Hamaker surfaces, and generally melt. In detail, they may exhibit 
either melting and nonmelting, depending on the metal and on the crystallo- 
graphy orientation. A minority of close packed surfaces such as Pb(lll) [3], 
Al(lll) and Al(100) [55], etc. display NM, and remain in fact smooth and dry 
all the way to T m . The vast majority of metal faces consisting of all other orien- 
tations, where packing is poorer, undergo SM. In intermediate packing cases 
like for instance Pb(100) [3], the first few layers melt, but the wetting does 
not proceed and the liquid film growth is blocked to a finite thickness until 
T m (incomplete SM). Although this incomplete melting will macroscopically 
appear just the same as NM — namely partial wetting with a finite wetting 
angle — there is a clear microscopic difference because the surface is no longer 
dry below T m . 

An early review of SM on metals can be found elsewhere [52]. We will limit 
ourselves to mention here a few high-temperature metal surface phenomena, 
particularly those connected with NM, which have been especially highlighted 
in our group. 

The first is the existence in surface NM of a critical liquid nucleation thick- 
ness £ cr u which is finite above T m and only vanishes at some "surface spinodal 
temperature" T s > T m . This fact, first discovered in simulation [2], indicates 
that in NM the solid surface is in fact metastable, and protected by a nu- 
cleation barrier, between T rn and T s . An immediate consequence is that NM 
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solid surfaces can be overheated, of course only for a short time, up to ati 
most T s > T m . For the (111) surfaces of Pb, Al and Au the calculated amount 
of maximum theoretical surface overheating T s — T m is 120 K, 150 K, and ~ 
150 K respectively. This possibility was demonstrated experimentally by laser 
heating techniques [56], and also in small Pb clusters [57]. 

Another phenomenon discovered theoretically is non-melting induced surface 
faceting. Consider a general crystal surface, whose orientation is close to, but 
not exactly coincident with, a flat NM face. Such a "vicinal surface" will consist 
of a sequence of flat terraces separated by surface steps. At low temperature 
the steps will usually repel mutually and form an ordered array, so that the 
solid vicinal surface is stable. As T m is approached, the steps suddenly coalesce 
in bunches giving way to much larger flat terraces, separated by very inclined 
facets where the steps have bunched up. The inclined facets with high step 
density actually melt, while the step-free flat facets remain dry. This non- 
melting induced surface faceting was first predicted thermodynamically [58], 
as seen on the free energy sketch (Fig. 14), demonstrated by MD simulation 
(Fig. 15) [59] , and eventually observed experimentally (Fig. 16) [60] . It is most 
likely the reason why e.g., beautifully perfect flat faces can be generated in 
gold, by simple heating close to the melting point. Au(lll) being a nonmelting 
face, it will at high temperature cause all steps and imperfections to bunch 
up in some place, sweeping itself clean and flat. 

A third interesting point is the microscopic investigation and clarification of 
the relationship between solid surface NM and its partial wetting by a drop of 
melt. While macroscopic partial wetting with a nonzero contact angle of the 
liquid with its own solid is not commonly reported for metals, it is very clearly 
found by simulation. Fig. 17 shows the contrasting fate of a simulated Al liquid 
nanodroplet at T = T m when approaching alternatively a SM Al(110) surface, 
or a NM Al(lll) surface. In the first case, the droplet spreads rapidly and 
disappears. In the second it settles down to a metastable state, forming the 
nanoscopic equivalent of a partial wetting angle 6lv = (22 ± 3)°. 

There is in fact, in a simple model, a direct relationship between 9 LV and the 
overheating O = T s — T m of a NM surface, namely 



where L is the latent heat of melting per unit mass, pi is the density of the 
liquid, £ is the correlation length in the liquid. The background of this formula 
is Eq. (2) with a short range interface interaction 



T 

— = 1 + 
T 




(7) 



V(£) = V(0)exp-(£/O 



(8) 
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plus the appropriate Young equation [5,58] 



Isv = Ilv cos 6 LV + 7sl cos 6 SL (9) 



where moreover 8sl ~ at T = T m . Based on Eq. (7) the calculated ideal 
overheating temperatures of Pb(lll), Al(lll) and Au(lll) given earlier above 
predict partial self-wetting angles 9lv of (16 ± 1)°, (18 ± 2)°, (33 ± 2)° respec- 
tively. By comparison Q^y is found to be (14 ± 1)° [60] in the experiment for 
Pb(lll), and (22 ± 2)° in the simulation for Al(lll) (see above). 

In addition to the above thermodynamics, there is a simple microscopic un- 
derstanding of the physics that leads to NM of these close packed metal sur- 
faces [5,31]. It has to do with the intrinsic atomic structure of the two SL and 
LV interfaces involved. A series of z-resolved in-plane averaged density pro- 
files of the SL and LV interfaces as displayed by MD simulations of suitably 
prepared liquid films on Au(lll), Au(110), and LJ(lll) near their respective 
T = T m is shown in Fig. 6. A well defined damped density oscillation is seen 
to propagate from the solid into the liquid, carrying approximately the solid 
interplanar distance as a wavelength. This wavelength of course depends on 
the crystallographic direction; in the present example it is large for (111) but 
small for (110). A second damped density oscillation starts at the liquid sur- 
face, carrying inward a generally different wavelength, determined this time 
by the main peak in the liquid structure factor S(k). This second oscillation, 
essentially non-existent in the LJ liquid, is generally strong in a metal (see 
discussion further below). Being a property of the liquid, the surface layering 
oscillation of the surface melted film is obviously face independent. 

When, as in Au(110), the solid and liquid oscillations facing one another pos- 
sess wavelengths that are out of tune, then their superposition is unfavorable 
and causes interface repulsion, eventually leading to SM. This is what hap- 
pens for Au(110). When instead, as in Au(lll), the two oscillations are close 
to being perfectly tuned, their superposition is favorable and causes interface 
attraction, leading to NM. Finally, the practical absence of a layering oscilla- 
tion in the LJ liquid surface indicates indifference of the two interfaces. They 
nonetheless eventually interact via the positive Hamaker constant [Eq. (5)], 
and therefore SM ensues in that case too. 
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4 Special cases 



4-1 Strained solids 

In an infinite bulk solid, a shift AP of external pressure provokes a shift of 
the bulk melting temperature by an amount AP(vi — v s )/(si — s s ), where vi 
and v s are the molar volumes of the liquid and the solid phase respectively, 
and si and s s are the corresponding molar entropies. Thus an increase of the 
external pressure generally increases the melting temperature. Contrary to 
this, a uniaxial strain always works the other way around and favors melting, 
irrespective of the sign of strain. Uniaxial compression or stretching increases 
the elastic energy of the solid (so long as the solid can sustain the correspond- 
ing stress without relevant plastic deformation). On the other hand the liquid 
does not support shear, and all elastic energy is released by flow. The bulk 
melting temperature of a strained solid is thus lowered by an amount which 
is quadratic in the anisotropic strain. This effect has been experimentally ob- 
served in He crystals [61]. In metals it has been studied theoretically through 
simple thermodynamics, and accessed directly by molecular dynamics simula- 
tions for aluminum [62,63]. This is therefore really a bulk effect though it will 
readily show up at surfaces, as shown by the simulation results of Fig. 18. Due 
to the strain-induced lowering of the bulk T m the liquid film thickness at fixed 
temperature on Al(llO) increases strongly and quadratically with strain, due 
simply to strain induced lowering of T m . A corollary of some practical impor- 
tance is that in any study of surface melting the strain conditions must be 
severely controlled. The Al test case indicates that a relatively small strain of 
10~ 3 implies a shift of T m of 0.03 degrees, or 0.003%. As such a degree of preci- 
sion is sometimes approached in surface melting studies, the error introduced 
by strains can be significant. Since strain conditions are seldom specified ex- 
plicitly, it seems possible and in some cases likely that some of the asymptotic 
SM data in the literature might need re-checking against that source of error. 

Besides this bulk effect, there are more interesting strain- related surface ef- 
fects. The first is the so-called Asaro-Tiller-Grinfeld instability, an experimen- 
tally well established effect; the second, is a strain-induced prewetting and 
eventual NM of an initially SM surface, an effect so far only hypothetical, but 
observed in MD simulation. 

The Asaro-Tiller-Grinfeld instability [64,65] is a periodic wavy modulation of 
the interface between the liquid and the strained solid, once the in-plane strain 
exceeds a critical threshold. It is connected with the possibility, permitted by 
a wavy interface only, of the solid "promontories" sticking out into the liquid 
to relax, reducing their strain and the associated elastic energy. When the 
bulk strain magnitude is large, the local surface strain release becomes large 
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enough, and pays for the free energy cost associated with the extra interface 
area due to the waviness. This effect has been observed at the He solid-liquid 
interface by the Ecole Normale group [61]. 

Strain induced prewetting can instead occur because a SV interface will gen- 
erally possess, in full equilibrium, a nonvanishing surface stress. For instance 
most metal surfaces are known to possess a large tensile surface stress, signi- 
fying that the topmost layers would really like to contract relative the bulk, 
whose lattice spacing is of course fixed [66] . When an in-plane bulk strain is 
applied, surface stress will work against the strain, or along with the strain, 
and the work reversibly done or released will increase or decrease respectively 
the overall SV surface free energy "fsv- As seen in Eq(l) and (2) a strain- 
induced increase of jsv wm further encourage SM, but a decrease will oppose 
it. In the latter case it can be shown (Fig. 19) in the very simple model of 
ref. [62] that strain will cause ordinary, continuous SM to be replaced by a 
prewetting transition: as T grows surface melting ceases to be continuous, and 
develops a sudden jump in the liquid film thickness I from zero to a finite value, 
from where it then grows continuously and diverges at T — > T m as in regular 
SM. This strain-induced prewetting transition has recently been confirmed by 
simulations [63], but has yet to be experimentally pursued. 



4-2 Defreezing of glass surfaces 

We dealt thus far with surfaces that are not just solid, but also crystalline. 
Many solids are however not crystalline, but amorphous or glassy. A glass 
does not break translational invariance the way a crystal does. Still, it breaks 
ergodicity in exactly the same manner as the crystal. Due to freezing and 
to the ensuing dynamical arrest, at the bulk glass temperature T g a glass, 
similar to a crystal, ceases to explore the full configurational phase space in 
the extensive way a liquid does, and acquires rigidity [67], the hallmark of 
broken positional ergodicity. 

Let us consider a semi-infinite glass with a free surface, and let T increase: 
will the glass surface defreeze at the same temperature T g as the bulk glass? 
In principle there is no simple answer to that, because the transition to the 
glassy state is a matter of dynamics, and not of thermodynamics. One can 
nonetheless build a simple "enthalpy functional approach", different in spirit 
but relatively parallel in reasoning, to the density functional theory of SM of 
Sec. 2. This approach is in turn based on a recent thermodynamic formulation 
of the glass transition, where this transition is in fact second order [68,69]. 

The enthalpy functional analog of Eq. 4 [70] which represents the surface ver- 
sion of the thermodynamical glass theory indeed predicts a surface defreezing 
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beginning at a lower temperature than the bulk T g , with a penetration depth 
that also grows logarithmically as T — > T g , superficially similar to SM in 
crystal surfaces. At a closer look, surface glass defreezing is in reality a much 
weaker affair than crystal SM is. To illustrate this, one can note a very substan- 
tial difference in the amount of total surface diffusivity, given by the integral 
of surface diffusivity over depth. This integral, measuring the amount of real 
flowing liquid at the surface, is divergent at T m for a SM crystal surface, but 
remains finite up to T g at a glass surface. 

Glasses may be expected to flow very, very slowly. If the flow takes place in 
the bulk, a piece of glass with a given initial shape will generally change its 
shape in a well defined manner [71]. The expected shape change will be totally 
different if the flow takes place at the surface instead of the bulk [72] . Unlike 
bulk flow, surface flow should wipe out sharp corners and edges, rounding 
off a finite but increasing length away from the corner as a function of time. 
It would be very interesting if these predictions, so far purely theoretical, 
could be explored experimentally, e.g., by studying shape changes of glasses 
by centrifuge experiments. 



4-3 Rotational surface melting in molecular solids 

We have been concerned with surfaces consisting implicitly of either atomic el- 
ements, or of small molecules such as NaCl merging together to form the solid. 
Let us consider by contrast molecular crystals made up of molecules that do 
not merge when crystallizing, and behave so to speak as nanosized rigid bodies, 
both in the crystal and in the liquid phase. In that case thermal disordering will 
concern, besides the positional degrees of freedom of the molecular centers of 
mass, also the rotational degrees of freedom. For very asymmetrical molecules, 
it can be expected that the two sets of degrees of freedom will be closely en- 
tangled. For nearly spherical molecules on the other hand the rotational and 
positional phenomena will involve very different energies and therefore differ- 
ent temperatures, the rotational well below the positional. In that situation 
one can have rotational melting, a first order phase transition transforming a 
proper crystal into a "plastic crystal" (somewhat of a misnomer) where the 
molecular centers of mass are still positionally crystalline, but where the rota- 
tional coordinates are thermally disordered, or melted. A particularly popular 
example of this behavior is provided by C 60 fullerite, where rotational melting 
is known to take place at T rot ~ 250 — 260 K [73] , whereas there is in Cqq no 
positional melting, and essentially no liquid phase. 

Consider the surface of such a molecular crystal and let T increase towards 
T rot . Will there be a rotational equivalent to surface melting, the first few layers 
thermally disordered and the bulk still rotationally ordered? The answer is in 
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principle yes; and again, simple free energy functionals can be resorted to in 
order to explore that scenario [74]. Here, perhaps the most interesting aspect 
here, (and certainly the most relevant to C 60 , where an abrupt surface phase 
transition was detected some 30 degrees below T rot [75]) appears to be the 
earliest stages of surface rotational disordering. 

There is a crucial difference between rotational and positional melting. Even 
in a rotationally melted plastic crystal, a molecule is still surrounded by a 
lattice, its cage exerting onto the molecule a sort of crystal field potential. In 
the ordinary liquid instead there is no remaining lattice at all. In the molec- 
ular solid, the crystal field potential felt by a first layer surface molecule is 
very different from the bulk one. The early stages of surface disordering and 
their precise nature will be vastly determined by that difference. Calculations 
show that at the surface of Cqo the crystal field actually frustrates bulk or- 
der; moreover it does so differently for different surface molecules. This gives 
in fact rise to a first layer thermal disordering process that proceeds in two 
stages, the more frustrated molecules disordering at a lower temperature, the 
less frustrated at higher temperature (but still below the bulk T rot ) [76]. 

Similar scenarios could be expected to take place for other molecular crystals, 
where surface studies are however not yet abundant so far. 



5 Liquid metal surfaces 

Our focus throughout has been on solid surfaces, touching only incidentally 
upon the physics of liquid surfaces, whenever that was important for surface 
melting. Yet, liquid surfaces have so to speak a life of their own. Liquid metal 
surfaces were long known to be interesting and structured. As revealed by e.g., 
X-ray reflectivity experiments [77,78], and as discussed long ago by chemical 
physicists such as Rice [79], the density profile of liquid metal surfaces shows a 
tendency to layering. Figs. 20 and 21 show as an example the density profiles 
obtained by simulation for the Pb, Al, LJ and NaCl liquid surfaces at their 
respective T m . 

Layering is in principle easy to understand, and simply represents the density 
response of the liquid to vacuum. Like all response functions, it is intimately 
related to the internal structure of the liquid, very thoroughly described in 
standard textbooks [80]. In short, if g(r) is the liquid pair correlation function, 
it will exhibit a main (damped) oscillation with some typical liquid interpar- 
ticle spacing A. Its Fourier transform is the liquid static structure factor S(k), 
which is correspondingly peaked around ko = 2tt/X. When an external pertur- 
bation, say a 5-function point-like repulsion, is inserted, the liquid density will 
locally drop, and will recover away from the point with a damped oscillation 
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dictated by S(k), and with a wavelength close to A. The surface represents 
another perturbation of that type (now extended) for the liquid metal, and 
layering is the result. Of course, the perturbation represented by vacuum is 
small only deep below the liquid surface, but is very large in the outermost 
layer. Here the precise form of layering will generally take a different form 
from that suggested by perturbative considerations. 

The metal atoms in the outermost surface layer are roughly speaking confined 
in two dimensions. Moreover, in some cases — like in the heavy noble metals 
that are inclined to surface reconstructions — the packing of surface atoms will 
tend to be closer and tighter than the corresponding bulk atoms. The question 
therefore arises whether the outer surface of a liquid metal could in some cases 
go as far as crystallizing in two dimensions (2D). Such a crystallization would 
represent a case of "surface freezing" — the opposite of surface melting — a phe- 
nomenon demonstrated and characterized in alkanes [81,82]. The additional 
interest of pursuing this in elemental liquid metals would be the possibility 
to use this for a particularly delicate test of the theories of two-dimensional 
freezing [83]. 

In two dimensions, freezing takes place in two steps. Disclination pairs should 
bind into dislocations with a first transition, implying the onset of a "hexatic 
phase" with power law orientational order. Upon further cooling, a second 
transition binds dislocation pairs, with the onset of power law positional or- 
der. These phenomena have been well documented in the melting of a 2D 
colloidal crystal [84]. A liquid metal surface might be able to exhibit, either 
above T m or in the supercooled state, a hexatic phase. Floating on top of a 
three dimensional liquid, the layered, nearly 2D solid surface would be able 
to exchange freely atoms with it, and disclinations/dislocations could readily 
form and dissolve avoiding all the delicate kinetic problems presented by strict 
2D systems. 

Celestini et al [85] carried out extensive MD simulations of the liquid Au sur- 
face as a function of temperature particularly in the supercooling regime. As 
Fig. 22 shows, an increasingly strong layering is indeed expected, and the top 
layer strongly resembles a 2D close packed lattice with a large number of discli- 
nations. Upon supercooling, disclinations rarefy and the system approaches a 
hexatic transition. On the brink of that transition however, bulk crystallization 
suddenly takes place, starting precisely from the liquid surface inwards [86]. 

An experimental attempt at pursuing experimentally the possible formation 
of 2D hexatics at liquid metal surfaces would seem very desirable, particularly 
with powerful tools such as surface X-ray reflectivity. Unfortunately, it appears 
that high melting point metals such as Au, Pt and Ir, where layering and 2D 
crystallization tendencies are particularly strong, are technically out of reach 
at least with this technique [77] . 
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6 Premelting of nanoclusters and nanowires 



Nanosystems became only recently a hot subject, but their thermal behavior 
has long been an important one. Small clusters liquefy readily or "premelt" , 
well below the bulk melting temperature of the same material. Basic thermo- 
dynamics predicts phenomenologically that a cluster of radius R should melt 
at T m (R) 

T m {R) « T m {l - -^r-[isv - (iSL + ILV)}} (10) 



This behavior was well verified experimentally in Au by classic cluster beam 
experiments by Buffat and Borel[87], whose results are shown in Fig. 23. 
A very parallel melting behavior is expected for one-dimensional extended 
nanosystems, such as nanowires [88]. In that case that Eq. (10) remains valid, 
although without the factor 2, as has also been confirmed by simulations. 

Microscopic theories as well as atomistic MD simulations of cluster melting are 
numerous and well documented. Reviews are given in standard books [89,90]. 
For Pb clusters, the size-induced drop of T m (R) is clearly recovered by simu- 
lations (Fig. 24). For Au clusters, of special experimental interest, MD simu- 
lations were carried out in order to understand how the premelted state was 
reached for increasing temperature. It was found first of all that the experi- 
mental size-driven lowering of T m (R) for Au clusters is confirmed. It can be 
seen in Figs. 25 and 26, that liquid-like diffusivity sets in at the surface below 
the cluster melting temperature (unless the cluster is so small as to melt at 
T m (oo)/2 or below). Melting at T m (R) occurs by a sudden propagation of this 
liquid "skin" into the bulk-like interior. This behavior is quite similar to that 
described by a phenomenological model of Celestini et al [91]. 

One would expect that clusters entirely bounded by NM crystal faces should 
not premelt. While some overheating has been confirmed in Pb clusters bounded 
entirely by (111) facets [57], this is not generally true. For example, clusters of 
systems with H < such as In [94] have been shown to premelt. The problem 
is that Eq. (10) is macroscopic, and as such it ignores microscopic details, such 
as edges and corners that clusters nonetheless possess. 

The simulated high temperature behavior of a NaCl cubic nanocluster ( a 
"nanogram of salt") is quite revealing [96]. The (100) nanocube faces are NM, 
and if they were unbounded, there would indeed be no premelting. However, 
the cube corners are the weak spot, and begin to melt even below T m . NaCl 
corner roughening was described long ago [95] however in a rigid step model, 
that did not permit melting. It remains in principle an open question how 
the picture should evolve with increasing grain size. Clearly the corners will 
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represent in all cases point-like germs of the liquid which are present below 
T m . While these germs will certainly enlarge into liquid pools, it remains to be 
clarified how the pools will spread out, and if and in what form their presence 
will permit, the faces of a large NaCl cube to exhibit NM and to remain 
solid up to and above T m . More generally, this approach should probably shed 
light also on the premelting of nanoparticles made of materials with H < 
mentioned above. 



7 Summary 

We reviewed the high temperature behavior of solid surfaces, in particular 
their tendency to wet spontaneously close to the bulk melting point. The case 
for which this does not happen ( and the melt somewhat unexpectedly fails to 
wet its own solid — surface nonmelting) has been treated in considerable de- 
tail. The large variety of behavior shown by different classes of solids has been 
addressed. We have also discussed other surface melting phenomena includ- 
ing rotational melting, glass surface defreezing, premelting of nanosystems. 
Finally, the existence of order at liquid surfaces has been highlighted as an 
important theme that is attracting increasing interest. It is hoped that this 
panoramic overview, despite its deliberate conciseness, will be of interest not 
only to the physics community, where it originates, but also to the materials 
science and chemical physics community, where it might be of some use. 
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Fig. 1. Melting of a pure and homogeneous material by heating along the solid- vapor 
coexistence curve up to the triple point. 
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Fig. 2. Experimental evidence of surface melting of Pb(llO): Medium energy ion 
scattering shows an increasing number of disordered (liquid) layers approaching the 
melting temperature. From Ref. [1]. 
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Fig. 3. Partial wetting of a solid by its own melt is an indirect evidence of surface 
nonmelting. 
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Fig. 4. (a) Minimizing grand potential (Eq. 4), (b) Mechanical equivalent (see text), 
(c) Density profile, showing that the SV interface has split into a SL plus LV inter- 
faces. 
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Fig. 5. Surface melting of LJ(llO): density profile and cristallinity profile. From 
Ref. [19]. 
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Fig. 6. Typical solid-liquid- vapor density profiles obtained by simulations of 
Lennard-Jones(llO), (111) and Au(llO) (three surfaces that undergo SM) just below 
T m . For Au(lll) (which is instead a NM surface) a nonequilibrium configuration 
is shown with two facing equal-period density oscillations which lead to attraction 
and eventually to the collapse of the two interfaces (see text). From Ref. [31]. 
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Fig. 7. Reentrant layering during the deposition of Argon on a graphite substrate. 
On the increasing T, the ellipsometric coverage isoterms switch from a staircase 
profile (growth layer-by-layer) to a continous one at 68. 8K, and then again to a 
staircase profile as temperature is further increased. The new structure consists of 
half layers, instead of full layers. Image and data from [41]. 
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Fig. 8. Top view (snapshot) of the three outer layers of simulated grand canonical 
Ar(lll) at 0.8 T m . This surface is nearly "disordered fiat". Atoms are represented 
as black, grey or white if their z coordinate indicates that they belong respectively 
to the subsurface layer, the surface layer or the adatom layer. For convenience, four 
adjacent simulation cells are shown. The nearly half occupancy of the surface layer 
(layer 1) is realized through large islands and large craters. Image from Ref. [46]. 



35 



Fig. 9. Argon bubble studies of liquid NaCl in contact with the solid. The 
solid-liquid- vapor junctions A and B reveal a surprising lack of complete wetting, 
with a large partial wetting angle of about 48 degrees. 
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Fig. 10. The simulated time evolution of NaCl liquid nanodroplet brought into the 
contact with NaCl(lOO) at the melting point. After 100 ps the drop stabilizes in 
the metastable state forming a partial wetting contact angle 6 = (50 ±5)° From 
Ref. [50]. 
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Fig. 11. Density profile of the liquid surface of NaCl obtained from a simulation at 
T m . The density drops very smoothly from the bulk in the liquid phase to that of 
the vapor one. From Ref. [96]. 
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Fig. 12. Schematic Ga phase diagram, showing the semimetallic a-phase surrounded 
by fully metallic liquid and solid-II phases. The picture for valence semiconductors 
such as Si, Ge is similar. 
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Fig. 13. (x, y) averaged atom density profile for (a) the Ge(lll)c(2x8) reconstructed 
surface at T = OK, (b) the same Ge(lll) surface at T ~ Tm, displaying "blocked 
melting" restricted to first bilayer and (c) crystalline bulk Ge at T ~ Tm- Note 
in (b) the disordering of the first bilayer with atoms distributed between different 
sublattices (inset) From Ref. [30]. 
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Fig. 14. Plot of a system with a solid (non-melting) branch and a liquid one. 
p = tan 9 and / = 7/ cos 9, where 9 is the tilting angle. The orientations up to pd 
remain dry, in the range (pd,Pm) there is non-melting induced faceting, and above 
p m surface melting occurs, pd might coincide with 0. From Ref. [58]. 
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Fig. 15. Non-melting induced faceting trajectories of the atoms in MD simulations, 
projected on a plane, (a) At low temperatures the Pb(432) surface is flat and stable, 
(b) At 0.97 T m the surface is phase separated in dry and wet facets. The atoms in 
the liquid part are recognizable by their wandering trajectories. From Ref. [59]. 
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Fig. 16. Non-melting induced faceting of Pb(lll) vicinals from Ref. [60]. Experi- 
mentally determined orientation angles as a function of T. Triangles refer to a tilt 
toward the (110) orientation, while circles refer to a tilt toward the (100) orienta- 
tion. Open symbols refer to the angle below which the vicinal is dry, while filled 
symbols to the angle above which the surface is melted. We refer to the original 
work for an explanation of the fit (solid and dashed curved). 
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Fig. 17. Evolution of Al liquid drop on a surface of bulk Al. Left column: drop on 
a surface undergoing surface melting (Al(110) at T = 0.99 T m ). (a) before contact; 
(b) after contact, the drop spreads readily; (c) the drop has been fully absorbed. 
Right column: drop on a nonmelting overheated surface (Al(lll) at T — 1.01 T m ). 
(d) before contact; (e) after contact: the drop settles but does not spread; (f) final 
drop shape. Darkness of atoms proportional to their square displacement in the run. 
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Fig. 18. Single snapshots (lateral view) of three simulated Al(llO) surfaces at the 
same temperature but under different strain conditions. Samples' size: 14 x 20 x 16. 
The atom stacking in the 20-atom rows orthogonal to the picture clearly distin- 
guishes the liquid and the solid phases. 
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Fig. 19. Phase diagram of a metal surface in presence of in-plane strain. Note that 
the negative (compressive) strain can give rise to a prewetting transition. From 
Ref. [62]. 
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Fig. 20. Density profile of the liquid surface of Pb and Al obtained from a simulation 
at T m . The profile shows layering oscillations at the liquid- vapor interface. 
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Fig. 21. Averaged density profile of the liquid surface of LJ (Ar) and NaCl ob- 
tained from a simulation at T m . The density drops smoothly from the bulk value 
in the liquid phase to that of the vapor phase, and the layering oscillations are 
undetectable. 
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Fig. 22. Nearly hexatic top layer of simulated liquid Au both above and below 
(undercooling) the melting point T m = 1335K. Fivefold (green) and sevenfold (red) 
disclinations are pinpointed. The corresponding density profiles along the surface 
normal are also shown. The first peak corresponds to the outmost layer of atoms in 
the maps. From Ref. [86]. 
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Fig. 23. Experimental values of the melting point temperature of gold particles. 
From Ref. [87]. 
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Fig. 24. Melting temperature as a function of iV -1 / 3 (where N is the number of 
atoms) of simulated Pb clusters. The points correspond to the discontinuities in the 
caloric curves in the inset. The solid line is the linear fit to the points. The triangle 
corresponds to bulk melting. From Ref. [92] 
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Fig. 25. Au§79: T dependence of the diffusion coefficient for different shells, com- 
pared with that for the supercooled bulk liquid. Inset: Radial density distribution 
at T = 700, 972 and 1049K. 
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Fig. 26. Logarithm of the effective Debye- Waller factor as a function of T for AU219, 
A11477, and AU879, and T dependence of S(k), with k along the (110) direction, for 
Au 879 . From Ref. [93]. 
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